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Abstract 

Using the general factorization approach, we present a detailed investigation for the branching 
ratios, CP asymmetries and longitudinal polarization fractions in all charmless hadronic B — > 
VV decays (except for the pure annihilation processes) within the most general two-Higgs- 
doublet model with spontaneous CP violation. It is seen that such a new physics model only 
has very small contributions to the branching ratios and longitudinal polarization fractions. 
However, as the model has rich CP-violating sources, it can lead to significant effects on the 
CP asymmetries, especially on those of penguin-dominated decay modes, which provides good 
signals for probing new physics beyond the SM in the future B-physics experiments. 
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I. INTRODUCTION 



During the recent years, tremendous progress in B physics has been made through 
the fruitful interplay between theory and experiment. The precise measurements of the 
i?-meson decays can provide an insight into very high energy scales via the indirect loop 
effects of new physics beyond the standard model (SM), which makes the study of exclusive 
non-leptonic 5-meson decays of great interest. 

In the SM, the phenomenon of CP violation can be accommodated in an efficient way 
through a complex phase entering the quark-mixing matrix, which governs the strength 
of the charged- current interactions of the quarks. This Kobayashi-Maskawa (KM) [H 
mechanism of CP violation is the subject of detailed investigation in these few decades. 
However, its origin remains unknown as it is put into the standard model through the 
complex Yukawa couplings. Moreover, the baryon asymmetry of the universe requires new 
sources of CP violation. Many possible extensions of the SM in the Higgs sector have been 
proposed [3| , and it was suggested that CP symmetry may break down spontaneously . 
A consistent and simple model, which provides a spontaneous CP violation mechanism, 
has been constructed completely in a general two-Higgs-doublet model (2HDM) 0, 
without imposing the ad hoc discrete symmetry, which is now commonly called as type 
III 2HDM. The type III 2HDM, which allows flavor-changing neutral currents (FCNCs) at 
tree level but suppressed by approximate U{1) flavor symmetry, has attracted much more 
interests. It is known that FCNCs are suppressed in low-energy experiments, especially 
for the lighter two generation quarks. Thus, the type III 2HDM can be parameterized in a 
way to satisfy the current experimental constraints. On the other hand, constrains on the 
general 2IIDM from the neutral mesons mixing (K^ — K^, — D^, and — B^) [g], 0] 
and from the radiative decays of bottom quark [8| have also been studied in detail. 

In recent years, there are many works about the i?-meson decays within the two-Higgs- 
doublet model. In Refs. jo], 10], the authors have studied the B PP,PV decays (with 
P and V denoting the pseudoscalar and vector mesons, respectively) within the type HI 
2HDM. Since through the measurements of magnitudes and phases of various helicity 
amplitudes, the charmless hadronic B VV decay modes can reveal more dynamics 
of exclusive B decays than B PP and B PV decays, in the present work we 
are going to make a detailed study for B VV decays within the type HI 2HDM by 
emphasizing on the new physics contributions. It will be seen that this speciflc new 
physics has remarkable effects on CP asymmetries, especially on the parameter Sf for 
the penguin-dominated decay modes. On the other hand, the new physics is found to 
have very small contributions to the branching ratios and the transverse polarizations. 
Furthermore, the polarization anomaly observed in B ^ pK* and B — ^ (pK* modes can 
not be improved in our current considered parameter spaces. 

The paper is organized as follows: In section II, we flrst describe the theoretical frame- 
work, including a brief introduction for the two-Higgs-doublet model with spontaneous 
CP violation, the effective Hamiltonian, as well as the decay amplitudes and CP violation 
formulas, which are the basic tools to estimate the branching ratios and CP asymmetries 
of i?-meson decays. In section HI, we list the Wilson coefficients and the other relevant 
input parameters. Our numerical predictions for the branching ratios, CP asymmetries 
and longitudinal polarization fractions are presented in Section IV. Our conclusions are 
presented in the last section. 
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II. THEORETICAL FRAMEWORK 



A. Outline of the Two-Higgs-doublet Model 

Motivated solely from the origin of CP violation, a general two-Higgs-doublet model 
with spontaneous CP violation (type III 2HDM) has been shown to provide one of the 
simplest and attractive models in understanding the origin and mechanism of CP violation 
at the weak scale. In such a model, there exists more physical neutral and charged Higgs 
bosons and rich CP violating sources from a single CP phase of the vacuum. These 
new sources of CP violation can lead to some significant phenomenological effects, which 
are promising to be tested by the future B factory and the LHCb experiments. In this 
paper, we shall focus on the phenomenological applications of the type III 2HDM on the 
two-body charmless hadronic B — > VV decays. 

The two complex Higgs doublets in the general 2HDM are generally expressed as 0, 



The corresponding Yukawa Lagrangian is given as 

'Cy = Vijai^i,L^aUj,R + ^ijai'i,L^aDj,R + h.C, (2) 

where the parameters rjija and ^ija are real, so that the lagrangian is CP invariant. After 
the symmetry is spontaneously broken down 

(0?)=t;ie^"S (0°) = ^;2e"^ (3) 

and the Goldstone particles have been eaten, the physical Higgs bosons are 

where are the charged scalar mass eigenstates, {(pi, ^s) are generally not the mass 
eigenstates but can be expressed as linear combinations of the mass eigenstates {H, h, A). 
Then the Yukawa part of the Lagrangian for physical particles can be written as 

Cy = rjli,i,LHiU,,R + rif^^Pi^LHiDj^n + i'^^i^i^AUj^R + ^giJiLH^Dj^R + h.c, (5) 

where 

Vij = Viji cos P + T]ij2e~^ sin P, 
^ij = -Viji(^~^ sin P + T]ij2 cos /3, 



'nfj = ^iji cos P + ^ij2e ^sinp, 

-D c -6. 



iij = -^ijie sinp + ^ij2Cosp, (6) 

and these couplings f]^ , tj^ , C,^ , are generally complex, which means CP violation. Ac- 
cording to the CKM mechanism, after diagonalizing the fermion terms' couplings rj^ and 
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T]^, the other couphngs become 

— TD^ — 771^ 

Ly = U^^—Un{v + c^l) + DL—Dn{v + (t>l) 

V V 

+Ud''DRH+ + DL~eDR{(t>l + i(t>l) + h.c, (7) 

i"" - vcKui''- (8) 

The Yukawa couphngs may be parameterized as foUowing 



with 



^ij — \j ■ (9) 



with V the vacuum expectation value v — 246 GeV. 



B. Effective Hamiltonian and decay amplitudes of B ^ VV decays 

Using the operator product expansion and the renormahzation group equation, the low 
energy effective Hamiltonian for charmless hadronic S-meson decays with AS = 1 can 
be written as 

= % E W.(Cl^3? + C2g^+ [C^Q^ + C[Q';\)+h.C., (10) 

p=u,c i=3,...,16 

where Cj(/i) {i = 1, . . . , 16) arc the Wilson coefficients that can be calculated by pertur- 
bative theory, and Qi are the quark and gluon effective operators, with Qi-io and Qn-ie 
coming from the SM and from the type III 2HDM, respectively. Their explicit forms are 
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defined as follows (taking b sqq transition as an example) 

Qi = {su)v-A{ub)v-A, 

Q2 = {SiUj)v-A{Ujbi)v-A, 

<33(5) = {sb)v-A^{qq)v-{+)A, 
q 

Qi(6) = isibj)v_A y^^{qjqi)v-(+)A, 
g 

^7(9) = 2(^^)^-^Xl^'?*^^^)^+(-)^' 
9 

Q8(io) = ■^{^ibj)v-A eg(g-,gi)y+(_)A, 

q 

Qii{i3) = {sb)s+py gg 5-(+)p, 

q 

r> '^qKq^\q) (- \ 

Wi2(i4) = \Sibj)s+p [qj qi)s-(+)p, 

q 

Qi5 = s-cr^^(l + 75)feV ^gcr^.(l + 75)g, 

q 

Q16 = + 75) bj J2 ^^^qj ^/..(l + 75) q^, (11) 

q 

where {qiq2)v±A = qilfi{l ± 75)'?2 and {qiq2)s±p = gi(l ± 75)'?2, with qu,d,s,c,b, and 
Cq is the electric charge number of q quark. The operators Qi in Eq. (fTTj) are obtained 
from via exchanging L ^ R, and we shall neglect their effects in our calculations 
for they are suppressed by a factor nig /nib in model III 2HDM. The Wilson coefficients 
Ci {i = 1, . . . , 10) have been calculated at leading order (LQ) [3, [ll] and at next-to- 
leading order (NLO) [l6| in the SM and also at LO in 2HDM [Ij, while Q (i = 11, . . . , 16) 



at LO can be found in Refs. 17, 18 



Having defined the effective Hamiltonian H^ff in terms of the four-quark operators 
Qi, we can then proceed to calculate the hadronic matrix elements with the generalized 
factorization assumption [l^, [13, 21, 2^ based on the naive factorization approach. 



For two-body charmless hadronic B — ^ VV decays, the decay amplitude of the local 
four fermion operators is defined as 

= ^{V,{h)V2{h2)\{q2q3)v±Aibqi)v^A\B), (12) 

where hi and h2 are the helicities of the final-state vector mesons Vi and V2 with four- 
momentum pi and P2, respectively. Since the B meson has spin zero, in the rest frame 
of i?-meson system, the two vector mesons have the same helicity due to helicity con- 
servation. Therefore three polarization states are possible in 5 — VV decays with one 
longitudinal (L) and two transverse, corresponding to helicities h = and h = ± (here 
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hi = h2 = h), respectively. We define the three hehcity amphtudes as follows 

A± = A{B ^Vi{puef)V2{p2,et)). 



(13) 



We choose the momentum p2 to be directed in the positive z-direction in the 5-meson 
rest frame, and the polarization four-vectors of the light vector mesons such that in a 
frame where both light mesons have large momentum along the 2;-axis. They are given 
by 



"1 

^1.2 



(0,±l,z,0)/v/2 



p^U/mi 



2, 



(14) 



where mi and m2 are the masses of Viand V2 mesons, respectively. Using the definitions 
for decay constants and form factors |23j, the tree- level hadronic matrix elements of the 
effective operators Qi can be decomposed as the following two amplitudes 



with 



(Fi(pi, e^,)\V - A\B){V2{P2, e^)|\/ - A|0), 
{Vi{pi,e\)\a^^{l + ^')\B){V2{p2,e'i)\a,.{l + ^'m. 



(15) 



(16) 



Here, for simplicity, we have omitted the quark spinors in the corresponding current 
operators in the above definitions. The three polarization amplitudes for and can 
be further written as 



Vi 



if ¥2^2 



A^'{mi+mB)TV^' 



2mB\pc\ 
rriB + mi 



0, 



2tf 



V2 



2T^''mB\pc\TT, 



^1 



m 



B 



mi 



(17) 



From the amplitude given by Eq. f|T5l) . the branching ratio for B —>■ VV decays then 
reads 

BriB VV) = {\Ao\' + \A^\' + \A^\') , (18) 



57rm 



B 



where tb is the lifetime of the B meson, and Pc is the center of mass momentum of either 
final-state meson with 



\Pc 



\/\m\ - (mi + m2)2] 



m 



B 



{rrii — m2Y 



2m 



(19) 



B 



In order to compare the relative size of the three different helicity amplitudes, we can 
define the longitudinal polarization fraction as 



h 





1^ 


2 


\Ao 


2 + 1 




2 + 


A_ 


2 



(20) 



which measures the relative strength of the longitudinally polarization amplitude in a 
given decay mode. 
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C. CP-violating asymmetries in i? — > VV decays 



Since there are abundant CP violation sources in the two-Higgs-doublet model, it is 
also necessary and interesting for us to discuss CP asymmetries in B ^ VV decays. 

Firstly, for charged i?^-meson decays, there is only one simple type of CP violating 
asymmetry, which detects direct CP violation 

T{B+^f+)-TiB'^f-) 



cv 



T{B+^f+) + TiB-^f- 



21 



For neutral S-meson decays, there is another type of CP violation coming from the 
mixing between B^ — B^ (here q = d or s) 



p 



g4t)\B') + g+lBl) 



[22) 



In this case, there are in general four amplitudes which can be expressed as 24, 25, 2^ 

(23) 



Af = {f\Heff\B',) , Af = {f\H^ff\B;), 
Af-={f\H,„\^^) , Af-={f\H,ff\B'^). 

For the B^ — Bd and Bg — Bg systems, the following approximations can be made 

(1 



both Bd and Bg systems : 



V 



only Bd system : AF ~ 0. 



(24) 



Using the decay amplitudes and the approximations listed in Eqs. (!23ll and (l24l) . the 
time-dependent decay probabilities for Bd system can then be written as 



nBl{t) f) : 

n'^M - /) 



\A,m+\\i 



\Afni + \\, 



e"^'{l + C/cos(Amt) - S'/sin(Amt)} , 
■e"^* {1 - Cf cos(Amt) + Sf sin(Amt)} , 



(25) 



while for Bg system, we have 



ns'git) - /) 



UfP(l + |A/ 



-Vt 



cosh 



/AFt 



D f sinh 



AFt 



-Cf cos(Amt) — Sf sin(Amt) 



cosh 



/AFt\ 



V 2 



Df sinh 



—Cf cos(Amt) + 5*/ sin(Amt) 



(26) 



where F is the average decay width, AF and Am are the width and mass difference, 
respectively. The other quantities are defined as 



A 



Cf = 



qAf_ 
pAf- 
1- |A/|^ 
1 + lAfP 



D 



f 



Sf 



2Re(Aj) 

i + iA/r 

21m(A/) 



1 + |A 



/I 



(27) 



From Eqs.( l25i) and ( 126|) . we can get: 



Acv{Bs 



f) = —Cf cos Amt + Sf sin Amt, 
—Cf cos Amt + Sf sin Amt 
^ cosh(^)+D^sinh(^^ 



(28) 



III. INPUT PARAMETERS 



The theoretical predictions in our calculations depend on many input parameters, such 
as the Wilson coefficients, the CKM matrix elements, the hadronic parameters, and so 
on. Here we present all the relevant input parameters as follows. 

It has been shown from B'^^ — ^ mixings that the parameters \\cc\ and \\ss\ in 
Eq. ffTTj) can reach to be around 100 |23], while their phases are not well constrained. In 
our present work we simply fix the phases to be 7r/4, and this choice will not cause any 
trouble in our numerical results. For the parameters At* and A^;,, the constraints come 
mainly from the experiments for B — B mixing, T{h — > S7), T{b — > ctVt), Pn, Rh, B — ^ PV 



13,118 



and the electric dipole moments (EDMS) of the electron and neutron [10|, 
Based on the above analyses, we choose the following three typical parameter spaces 
which are allowed by the present experiments and have been adopted for the B PV 
decays [F" 



Case 
Case 
Case 



A 
B 

C 





= 0.15; 


Abfe 


= 50, 


A«| 


= 0.3; 


Afcb 


= 30, 


A«| 


= 0.03; 


Afc?) 


= 100, 



and 6tt + 6bb = it/2 

responding to the SM, we use the results listed in the paper 
coefficients in the type III 2HDM, we redefine them as Cn 



For the Higgs masses and the Wilson coefficients of Ci,...,io cor- 

10 , while for the Wilson 



16 



Cii ... 16 in order to 



compare the contributions from those operators in SM, here the factor 



nit 



is associated 



with the operators in 2HDM, the numerical values for Cii,...,i6 are listed in Table [B 



TABLE I: The Wilson coefficients Cn, 
4.2 GeV in 2HDM. 



,16 



(*) 



mi, 



-c 



11,. ..,16 



in & — > s transition at fi 



rub 



Case A 



Case B 



Case C 



Cu 
C12 

Cl3 
Cl4 
Cl5 
C16 



-0.0085 + 0.012i 


-0.0030 - 0.0049i 
-0.000060 - O.OOOlOf 
0.000033 + 0.000055i 
-0.00010 - 0.00017i 



-0.0085 + 0.018i 


-0.0052 - 0.0069i 
-0.00011 - 0.00014? 
0.000058 + 0.000078i 
-0.00018 - 0.00024i 



-0.010 + 0.012i 


-0.0029 - 0.0052i 
-0.000059 - O.OOOlOi 
0.000032 + 0.000059i 

-0.0001 - 0.00018f 



As for the CKM matrix elements, we shall use the Wolfenstein parametrization [28 
with the values [26]: A = 0.8533 ± 0.0512, A = 0.2200 ± 0.0026, p = 0.20 ± 0.09, and 
f] = 0.33 ±0.05. 
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For the hadronic parameters, the decay constants, and the form factors, we hst them 
in Tables. Hi] and IIIH respectively. 



TABLE II: The hadronic input parameters 2g] and the decay constants taken from the QCD 
sum rules 29 ] and Lattice theory [sB] . 









Mb, 


Mb. 


nib 


1.638ps 


1.528ps 


L472ps 


5.28GeV 


5.37GeV 


4.2GeV 


mt 


rriu 


rrid 


nic 


rus 


rUpO 


174GeV 


3.2MeV 


6.4MeV 


l.lGeV 


0.105GeV 


0.77GeV 


mp± 


rriuj 




m^*± 


m^.o 


Aqcd 


0.77GeV 


0.782GeV 


1.02GeV 


0.892GeV 


0.896GeV 


225MeV 


fp 








Jp 


J UJ 


0.205GeV 


0.195GeV 


0.217GeV 


0.231GeV 


0.147GeV 


0.133GeV 


Ik* 


J(f> 










0.156GeV 


O.lSSGeV 











TABLE III: The relevant B ^ V transition form factors at = taken from the light-cone 



sum rules (LCSR) [3l|,|32| 



decay channel 


V 


^0 


Ai 


A2 


Ti 


Ts 


B^p 


0.323 


0.303 


0.242 


0.221 


0.267 


0.176 


B ^uj 


0.293 


0.281 


0.219 


0.198 


0.242 


0.155 


B ^ K* 


0.411 


0.374 


0.292 


0.259 


0.333 


0.202 


Bs^K* 


0.311 


0.363 


0.233 


0.181 


0.26 


0.136 


Bs^cp 


0.434 


0.474 


0.311 


0.234 


0.349 


0.175 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, we shall classify the 28 channels of , and Bg decays into two light 
vector mesons according to the reliability of the calculation for various observables, which 
is motivated by the dominated contributing operators. We shall give our predictions for 
the branching ratios, the CP asymmetries, and the longitudinal polarization fractions 
both in the SM and in the 2HDM. Comparisons with the current experiment data, if 
possible, are also made. 

Before moving to the detailed discussions, some general observations of new physics 
effects on B ^ VV decays should be made. As can be seen from Eqs. ([9]) and (ITTl) . 
the contributions of new physics operators Oii^...^i6 are always proportional to the factor 
nLq/v. Thus, they are severely suppressed for the first generation quarks. In this case, 
for B — > pK*,u!K*, pp, up, uu and Bg pK*,ujK*, K*K*, pcj), uxp decay channels, we can 
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safely ignore the contributions from those new operators. Note that the new physics still 
has effects on the Wilson Coefficients Ci_io. On the other hand, for B —>■ (j)K*,(j)p,(j)u) 
and Bg (j)K*,(f)(j) decay channels, since these are all induced by 6 — qss (g = d,s) 
transitions, we could not ignore the new operators' contributions any more in this case. 
In the general factorization approach, it is impossible to produce a vector meson via the 
scalar and/or pseudoscalar currents from the vacuum state, and hence the new operators 
Qii and Q13 have no contributions to -B — >■ VV decays. Moreover, from the results listed 
in Table [H it can be seen that all the contributing new operators Q12, 14,15,16 have only 
very small (even zero) Wilson coefficients. It is therefore expected that the new physics 
will have very small effects on the branching ratios and transverse amplitudes (hence on 
the transverse polarization fractions) of i? — > VV decays. 

A. CP-averaged branching ratios and direct CP violation 

According to different decay modes, we shall give our predictions for the branching 
ratios and direct CP violations one by one. 

(i), color-allowed tree-dominated decays. Our predictions for the CP-averaged branch- 
ing ratios and the direct CP asymmetries are presented in Table IIVI From the numerical 
results, we can see that the branching ratios are all at 10~^ order, and the direct CP asym- 
metries are all very small since the penguin amplitude contributions are much smaller than 
the ones from the tree diagrams. Most predictions within the SM are consistent with the 
current experiment data, and the new physics has very small effects on these types of 
decays. 

TABLE IV: The CP-averaged branching ratios (in unit of 10^^) (first line) and the direct CP 
violations (second line) for the color-allowed tree-dominant processes both in the SM and in the 
type III 2HDM. Case A-C stand for the three different parameter spaces listed in Section III. 



Decay modes Case A Case B Case C SM Exp. 



B+ ^ 


' P 


v° 


14.59 


14.59 


14.59 


15.53 


18.2±3.0 








-0.004 


-0.004 


-0.004 


-0.002 


-0.08±0.13 




P' 




26.33 
-0.043 


25.93 
-0.043 


26.73 
-0.042 


27.49 
-0.035 


94 0+3-1 
^^•^-3.2 


B+ ^ 


' P 




12.66 


12.47 


12.85 


13.97 










-0.042 


-0.043 


-0.042 


-0.034 


0.04±0.18 




P^ 


K*- 


36.88 
-0.043 


36.32 
-0.043 


37.44 
-0.042 


38.50 
-0.035 





(ii), color-suppressed tree-dominated decays. The numerical results are given in Ta- 
ble lYl it is interesting to note that the branching ratios will generally become smaller 
after including the new physics contributions except for the B — > p°p° mode. Further- 
more, there are big direct CP violations in these decay processes except for the B~^ — > p^u 
mode, and the new physics has more effects on the direct CP asymmetries than on the 
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branching ratios through the Wilson coefficient functions, although there are no new 
operator contributions to the hadronic matrix elements in this type decays within our 
approximations. Compared to Case A and Case C, Case B has the biggest corrections to 
the CP asymmetries of the SM. 



TABLE V: The same as Table ITVl but for color-suppressed tree-dominant processes. 



Decay modes 


Case A 


Case B 


Case C 


SM 


Exp. 






0.0814 


0.0897 


0.0754 


0.065 


0.86±0.28 






0.176 


0.218 


0.119 


0.153 




B+ - 


-> LOU) 


0.112 


0.110 


0.115 


0.160 


<4.0 






-0.117 


-0.088 


-0.144 


-0.207 




Bs- 




0.081 


0.090 


0.073 


0.092 


< 7.67 X lO-'' 






0.176 


0.218 


0.119 


0.153 




Bs- 




0.183 


0.180 


0.187 


0.262 








-0.167 


-0.088 


-0.144 


-0.207 




B+ - 


-^p^LJ 


0.024 


0.024 


0.024 


0.076 


<1.5 






-0.063 


-0.063 


-0.063 


-0.035 





(iii) , penguin-dominated decays. We may divide such decays into two types: AS* = 1 
and AD = 1 decay modes. They are corresponding to the upper and the lower parts 
in Table IVIt respectively. From the numerical results, we can see that all the eleven 
AS* = 1 decay modes have branching ratios up to 10~^ or even to 10~^ order, since 
they involve the relative large CKM matrix elements V^*, while the AD = 1 ones have 
much smaller branching ratios of order of 10"'' due to the smaller CKM matrix elements 

For B — >■ uK* and Bg — > 00 decay modes, our predictions for the branching ratios 
with including the new operator contributions have similar results as the ones within the 
SM, which, however, are not quite consistent with the current experimental data; the 
numerical results for B —>■ uK* modes are larger than the current experiment limit, and 
the prediction for Bs — *■ 00 is about two times larger than the present data. For the other 
decay modes, our predictions for the branching ratios are in general agreement with the 
data. As for the direct CP asymmetries, there are big CP violations in some decay modes, 
and the new physics can lead to remarkable effects. Our predictions are consistent to the 
data in all these decay modes. 

(iv) , electroweak penguin or QCD flavor singlet dominated decays. As can be seen 
from Table IVIH this type of decays are expected to have smaller branching ratios due to 
the large cancelations among the different Wilson coefficients. Although there are new 
operator contributions in B ^ p(f) and uj(f) decay modes, the predicted branching ratios are 
still small. The direct CP asymmetries for these decays are all small, and the new physics 
effects on these observables are not prominent. Due to the lack of accurate experimental 
data, we couldn't compare our predictions with the data yet. 

(v) , the pure annihilation decays. Only six decays belong to this class, namely B^ 
K*+K*-, B° 00, B, ^ p+p-, Bs -> pV°, p^u, and ^ uu. Due to the lack 
of the information for the Vi V2 transition form factor at large momentum transfers, 
we shall not consider them in details in this paper. 
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TABLE VI: The same as Table ITVl but for the penguin-dominated decay modes. The upper and 
the lower parts correspond to AS" = 1 and AD = 1 processes, respectively. 



Decay 


modes 


Case A 


Case B 


Case C 


SM 


Exp. 


B+ ^ 


■ p+K*'^ 


7.169 


7.409 


7.027 


7.287 


9.2±1.5 






0.084 


0.117 


0.049 


0.018 


-0.01±0.16 


B+ ^ 


■ p^K*+ 


5.853 


6.229 


5.526 


5.575 


<6.1 






0.184 


0.196 


0.169 


0.122 


n 90+0-32 






6.396 


6.513 


6.324 


6.245 


5.6±1.6 






0.054 


0.073 


0.033 


0.018 


0.09±0.19 




p-K*+ 


6.046 


6.738 


5.445 


5.571 


<12 






0.295 


0.301 


0.283 


0.199 






ujK*^ 


3.412 


3.513 


3.351 


3.498 


<2.7 






0.078 


0.107 


0.048 


0.024 






■ ujK*+ 


3.247 


3.5697 


2.965 


3.123 


<3.4 






0.265 


0.274 


0.251 


0.176 








9.276 


9.704 


9.221 


9.318 


9.5±0.8 






0.045 


0.081 


-0.002 


0.020 


-0.01zb0.06 




■ (j)K*+ 


9.867 


10.32 


9.775 


9.979 


lO.Oibl.l 






0.039 


0.074 


-0.013 


0.020 


-0.01±0.08 


Bs^ 


</,</. 


28.99 


30.34 


28.64 


28.85 


14^7 X 10-6 






0.054 


0.089 


0.006 


0.020 




Bs^ 




9.303 


9.614 


9.118 


9.456 


< 1.681 X 10-^ 






0.084 


0.117 


0.049 


0.018 




Bs^ 


K*+K*- 


8.404 


9.366 


7.569 


7.744 








0.295 


0.302 


0.283 


0.199 








0.410 


0.420 


0.413 


0.408 


0.49t°il 






-0.092 


-0.061 


-0.133 


-0.145 




B+^ 


■ K*^K*'^ 


0.439 


0.450 


0.443 


0.437 


< 2.2 






-0.092 


-0.061 


-0.133 


-0.145 




Bs^ 




0.517 


0.532 


0.521 


0.526 


< 1.013 X 10-3 






-0.094 


-0.056 


-0.145 


-0.161 





B. Time-dependent CP violating parameters Cf, Sf and Dj 

Since there are abundant CP violating sources in type III 2HDM, it is expected that 
there are relatively large CP violations in 2HDM than in the SM. Using the relevant 
formulas given in section II, we can predict the time-dependent CP asymmetries in neutral 
Bd and Bg decays, with the numerical results given in Tables IVIIII and IIXI respectively. 

From these two tables, it is seen that, for B^ p^p^ , p^(j) and a;0 decay modes, the new 
physics has hardly any effects on the parameters Cf and Sf, even though there are new 



12 



TABLE VII: The same as Talbe IIVI but for the electroweak pengum or QCD flavor singlet 
dominated decays. 



Decay modes 


Case A 


Case B 


Case C 


SM 


Exp. 


B+ p+(j) 


0.0054 


0.0054 


0.0054 


0.0043 


< 16 




-0.011 


-0.011 


-0.011 


-0.014 




B° p^cp 


0.0025 


0.0025 


0.0025 


0.0020 


<13 




-0.011 


-0.011 


-0.011 


-0.014 




B° ujcj) 


0.0022 


0.0022 


0.0022 


0.0017 


< 1.2 




-0.011 


-0.011 


-0.011 


-0.014 




Bs ^ 


0.796 


0.796 


0.796 


0.687 


< 6.17 X 10-^ 




0.0048 


0.0048 


0.0048 


0.0039 




Bs (puj 


0.038 


0.038 


0.038 


0.045 






0.020 


0.020 


0.020 


0.018 





operators contributions in B p (j) and ujcj) decay modes. On the other hand, the new 
physics has remarkable effects on the other decay modes, especially on uouj one (for 

this mode the new physics can even change the sign of the parameter Sf). Furthermore, 
different parameter spaces also have remarkable effects on these CP violation parameters. 

For Bs system, there are new operator contributions only in Bs (fxp mode. As is 
expected, the new physics has remarkable influence on the parameters Cf, Sf, and Df. 
For the other four decay modes, although there are no new operator contributions, the 
new physics still has big effects on the parameter 5*/, but small effects on C/ and Df. 



C. The polarization in B ^ pK* and (pK* decays 



Motivated by the polarization anomaly observed by the BarBar [33|, Belle [3J] and 
CDF js^ experiments, we shall study the polarization in B —>■ VV decays, especially in 
B — > pK* and (pK* decays in this section. 

One important point that should be noted is that the predictions for the branching 
ratios of -B ^ pK* and (pK* modes are well consistent with the experiment data, which 
means that if we want to solve the observed polarization anomaly, we need to find some 
way to reduce the longitudinal amplitude and enhance transverse ones simultaneously. 
Many studies have been made to try to provide possible resolutions to the anomaly both 
within the SM I 



371 . 138l . |39| and in various new physics models |40|, |41|, |42|. Here we 



only concentrate on the longitudinal polarization fraction and the main results are listed 
in Table IS 

It is noted that the polarization anomaly could be well resolved by introducing the 
tensor operators Oti = sa^'^{l+'y^)bsa^,y{l+'y5)s and Ors = Sia^'^ {l+'y^)b 



in Ref. [42j. It is interesting to see that these two operators have similar forms as Qis and 
Qi6 in Eq. (Ilip . However, from the numerical results given by Table IX] we can see that 
the predicted longitudinal polarization fraction fi for these decay modes in the type HI 
2HDM is almost the same as the one within the SM. Although there are new operator 
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TABLE VIII: The time-dependent CP asymmetry parameters Cf (first line) and Sf (second 
line) for decays both in the SM and in the type III 2HDM. Case A-C stand for the three 
different parameter spaces listed in Section III. 



Decay modes 


Case A 


Case B 


Case C 


SM 


B^ p+p- 


0.043 


0.043 


0.042 


0.035 




-0.95 


-0.95 


-0.95 


-0.95 




-0.18 


-0.22 


-0.12 


-0.15 




0.97 


0.92 


0.99 


0.89 


S° ^ up^ 


0.063 


0.063 


0.063 


0.029 




-0.61 


-0.61 


-0.62 


-0.97 


5° ^ 4>p^ 


0.011 


0.011 


0.011 


0.014 




0.70 


0.70 


0.70 


0.70 


B°^u(t> 


0.011 


0.011 


0.011 


0.014 




0.70 


0.70 


0.70 


0.70 


B^ ^uu 


0.12 


0.09 


0.14 


0.21 




0.53 


0.65 


0.40 


-0.18 




0.092 


0.061 


0.13 


0.15 




0.85 


0.92 


0.75 


0.57 



TABLE IX: The time-dependent CP asymmetry parameters Cf (first line), 5"^ (second line), 
and Df (third line) for Bg decays both in the SM and in the type III 2HDM 



Decay modes 


Case A 


Case B 


Case C 


SM 


Bs ^ 


-0.005 


-0.005 


-0.005 


-0.004 




0.052 


0.052 


0.052 


0.14 




0.99 


0.99 


0.99 


0.99 


Bs — > 4>UJ 


-0.020 


-0.020 


-0.020 


-0.018 




0.23 


0.23 


0.23 


0.49 




0.97 


0.97 


0.97 


0.87 


Bs^(i>ct> 


-0.054 


-0.090 


-0.060 


-0.020 




0.33 


0.49 


0.14 


-0.004 




0.94 


0.87 


0.99 


1.0 


Bs K*+K*- 


-0.30 


-0.30 


-0.28 


-0.20 




0.92 


0.95 


0.88 


0.79 




0.25 


0.12 


0.39 


0.57 


Bs K*^K*^ 


-0.085 


-0.12 


-0.049 


-0.018 




0.31 


0.45 


0.15 


-0.003 




0.95 


0.88 


0.99 


1.0 
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contributions in B ^ (pK* modes, we still can not resolve the polarization anomaly 
observed in this decay mode. This is due to the fact that the strength of new operators 
in 2HDM is severely suppressed by the factor mqXqq/rrib. Moreover, as has already been 
mentioned in the beginning of this section, the Wilson coefficients of these new operators 
are very small, which also result in the small effects on the transverse amplitudes. 

TABLE X: The longitudinal polarization fractions Jl for B pK* and (f>K* decay modes. Case 
A-C stand for the three different parameter spaces in the type III 2HDM. 



Decay modes 


SM 


Case A 


Case B 


Case C 


Exp. 


B+ p+K*° 


0.91 


0.91 


0.91 


0.91 


0.48 ± 0.08 


B° pK*° 


0.95 


0.95 


0.93 


0.93 


0.57 ±0.12 


B+ (pK* 


0.89 


0.89 


0.89 


0.89 


0.50 ± 0.05 


B° (t)K*^ 


0.89 


0.89 


0.89 


0.89 


0.491 ± 0.032 



For the other B — > VV decay modes, the predictions for longitudinal polarization 
fractions are always about 0.90 ~ 0.95. For simphfy, we shall not list the results in details 
anymore. 

In conclusion, adopting the current parameter spaces and with the general factorization 
method, we could not resolve the polarization anomaly observed in S — > pK* and (t)K* 
modes within the SM and 2HDM. 

V. CONCLUSIONS 

Using the general factorization approach, we have studied all the B VV decay 
modes except for pure annihilation decay channels both within the SM and in the two- 
Higgs-doublet model. From the numerical results given in the previous section, we can 
see that: for the branching ratios, our predictions are generally well consistent with the 
current experimental data expect for the B^ (fxf) decay mode, and the new physics has 
margin or even negligible effects on this observable. However, the new physics can give 
remarkable contributions to the CP asymmetry parameters C/ and Sf, especially to Sf in 
the penguin-dominated decay modes. Unfortunately, our predictions for the longitudinal 
polarization fractions of i? — > pK* and (j)K* decay modes in 2HDM arc still as large 
as the ones in the SM, which are much larger than the experimental data. Some new 
mechanisms may be needed to improve those discrepancies. 

For simplicity, in this paper we have neglected the contributions from annihilation and 
exchange diagrams, although they may play a significant rule in some decay channels. In 
our numerical calculations, we have only considered three possible parameter spaces for the 
type III 2HDM. Also we have totally neglected the first generation Yukawa couplings and 
the off-diagonal matrix elements of the Yukawa coupling matrix, in order to eliminate the 
FCNC at tree level. However, it is possible that the FCNC involving the third generation 
quarks still exists at tree level, making the constraints less stronger. In a word, we do 
not exclude the possibility to improve the predictions by using the other factorization 
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methods with the annihilation and exchange diagram contributions included, by choosing 
other parameters spaces, or even by introducing additional fourth-generation quarks [43| . 

In conclusion, we have shown that the new Higgs bosons in the type III 2IIDM with 
spontaneous CP violation can have significant effects on some charmless B VV decays, 
especially for the penguin-dominated decay modes, which can be used as good signals to 
test the SM and to explore new physics from more precise measurements in the future 
-B-factory experiments. 
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